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We have measured the complex conductivity, a, of a Bi2Sr2CaCu20s+s (BSCCO) thin film between 
0.2 and 0.8 THz. We find a in the superconducting state to be well described as the sum of 
contributions from quasiparticles, the condensate, and order parameter fiuctuations which draw 
30% of the spectral weight from the condensate. An analysis based on this decomposition yields a 
quasiparticle scattering rate on the order of ksT /% for temperatures below Tc. 



74.25.Gz, 74.25.Nf, 74.25.-q, 74.72.Hs 

The unusual properties of the quasiparticle (QP) life- 
time Tqp in high- Tc cuprate superconductors provide 
evidence for a non-Fermi liquid normal state and a 
non-BCS mechanism. Because angle-resolved photoe- 
mission (ARPES) shows that QP properties are highly 
anisotropic it is important to distinguish the be- 

havior of Tqp at different parts of the Fermi surface. Near 
the maximum of the d-wave gap at the (tt, 0) point in 
momentum space, the ARPES lineshape narrows rapidly 
below Tc . The change suggests the onset of a well- 
defined QP in this antinodal region. Recently, improve- 
ments in detectors have made it possible to resolve the 
lineshape in the nodal or (tt, tt) direction as well. Valla et 
al. Q reported that in the normal state, Tqp can be de- 
scribed by the marginal Fermi liquid phenomenology , 
l/r^p ~ max(tj,T). Surprisingly, l/r^p of nodal quasi- 
particles (at the Fermi surface) remains approximately 
proportional to T in the superconducting state as well 
Thus the change in the lifetime of such QP's upon 
entering the superconducting state is much less than that 
of the antinodal QPs. 

The nodal QP's apparent insensitivity to superconduc- 
tivity is surprising for two reasons. First, one might 
expect the reduction in the phase space for scattering 
which accompanies the formation of a gap to increase the 
QP lifetime. Second, Tqp, as determined by microwave 
[0-^ and thermal conductivity measurements in the 
closely related compound YBa2Cu30j-s (YBCO) does 
indeed increase rapidly below Tc. While such trans- 
port experiments are not momentum resolved, there is 
strong evidence that they probe the nodal QP lifetime. 
It was reported in Ref. that the microwave conduc- 
tivity, a — ai + ia2, can be described (at low T) by a 
two-fluid model which sums the contributions of conden- 
sate and QPs. The spectral weight in the normal fluid, 
Pn J fJidv, grows linearly with T at low temperature 
at a rate which matches the linear decrease in superfluid 
density ps{T), measured independently via the penetra- 



tion depth. The linear rate of transfer from ps to pn 
agrees with the theoretical prediction for nodal QPs in a 
d-wave superconductor pT[ |. 

Thus, the transport measurements of Tqp in YBCO 
stand in apparent contradiction to ARPES measure- 
ments in BSCCO. It would appear that either the single- 
particle lifetime probed by ARPES is qualitatively dif- 
ferent than the transport lifetime, or that Tqp differs in 
BSCCO and YBCO despite their similar bilayer struc- 
ture and Tc- The resolution of this issue requires inter- 
changing the two probes and the two material systems. 
Since ARPES lineshapes in YBCO are difficult to ob- 
tain, it is important to focus on the transport properties 
of BSCCO, for example the microwave conductivity 

Unfortunately, it has proved difficult to determine Tqp 
reliably from microwave data in BSCCO. In contrast to 
the strong dispersion in YBCO, cti is approximately con- 
stant in the range from i^=10-50 GHz [p"3|-p^. Conse- 
quently, the transport lifetime is not directly resolved 
in the conductivity spectra, although they do indicate 
that Tqp <C l/27r-50 GHz, which is much smaller than 
in YBCO. Despite the lack of frequency dependence, in 
principle Tqp could still be estimated from the magni- 
tude of ai, which (in the coTqp <^ 1 limit) is proportional 
to PnTqp. Using the two-fluid model, the magnitude of 
Pn(T) follows from measurement of ps{T). In practice, 
this analysis yields results which are not self-consistent, 
for the following reason. Just as in YBCO, the conden- 
sate density varies as Ps(0) — aT at low T p^-p^, im- 
plying that the samples are in the clean d-wave regime 
where pn = aT . However, ci in BSCCO does not tend 
to zero as T approaches zero. Even at the lowest temper- 
atures measured, about 5 K, ui remains approximately 
8 times larger than the normal state value just above Tc, 
'^n(Tc), and is far larger than the 'universal conductiv- 
ity' proposed by Lee [16|. Thus Tqp calculated from ai/pn 
appears to diverge, which is inconsistent with the obser- 
vation that the conductivity is frequency independent in 
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the microwave regime. To extract a finite Tqp one may 
assume that a large fraction of the quasiparticles remain 
uncondensed, i.e. Pn{T) = /O„(0) + aT [|l^, or subtract 
the zero temperature hmit of the dissipation prior to the 
two- fluid model analysis p^ ]. 

In view of the uncertainties inherent in these assump- 
tions, it is desirable to extend the conductivity mea- 
surements to the frequency range probed by ARPES, 
K, ksT/h. Here we report measurements of cr in a 
BSCCO thin film in the frequency range from 0.18 to 
1.0 THz, above the domain of microwave and below that 
of infrared spectroscopies. To cover this region, we used 
time-domain terahertz spectroscopy, a technique which 
is based on the generation and detection of single cy- 
cle electromagnetic pulses. The transmission coefficient 
of such pulses yields both the real and imaginary parts 
of cr directly and independently |l9| ]. The sample, with 
Tc=85K, was grown on a LaAlOs substrate using atomic 
layer-by-layer molecular beam epitaxy |20[ . The resis- 
tance R versus T is linear and R{0)/R{300K) = 4x 10-^ 
where R{0) is the extrapolation of the linear resistance 
to zero temperature. 

The terahertz data show a strong frequency depen- 
dence in a below Tc, from which Tqp in the supercon- 
ducting state can be determined directly. We find that 
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FIG. 1. Left panel: ai plotted versus T for 0.2, 0.3, 0.4, 
0.6 and 0.8 THz as squares, octagons, diamonds, circles and 
triangles. At each frequency the arrows mark the T where cti 
begins to decrease. Upper right panel: l/r^p implied by the 
arrows plotted versus T, with the same T scale as the panel 
below. The dashed line is a fit to l/r^p oc T. Lower right 
panel: Drude conductivity, aqp{T), using Tqp proportional to 
the linear fit in the upper panel plotted on the same scale as 
the left panel. The same frequencies are shown as the left 
panel, increasing v having smaller aqp{T). 

l/Tqp{T) « 0.8kB{T+ 10K)/n, consistent with ARPES 
measurements Q,^. Furthermore, we find the total con- 
ductivity to be incompatible with the two-fluid model. 
Instead, our measurements indicate that the conductiv- 
ity comprises three contributions: QPs, the condensate, 



and a low frequency collective mode with spectral weight 
drawn from the condensate. 

We now turn to ai(w,T) and discuss its salient fea- 
tures. The left panel of Figure 1 shows cti plotted versus 
T for a representative sample of frequencies in our range. 
At temperatures above Tc there is little frequency de- 
pendence, as we would expect since l/r^p ^ 2tt-1 THz 
in this temperature range. Below Tc, however, <ti{T) de- 
pends strongly on frequency. At i'—Q.2 THz the T depen- 
dence is quite similar to that measured in the microwave 
regime on BSCCO single crystals [p3|-p^. cti rises to a 
maximum at 25 K, but at 5 K remains much greater than 
Cn{Tc)- As the frequency increases from 0.2 to 0.8 THz 
the maximum shifts from 25 K to Tc, indicating that 
1/Tgp is sweeping through our bandwidth. For a given 
frequency, the QP conductivity will decrease below the 
temperature where l/rgp(T) k uj. Kt each measurement 
frequency, this T is marked with an arrow. The mea- 
surement frequency (in radians/s) and corresponding T 
for each arrow are shown in the upper right panel of Fig. 
1. This admittedly rough estimate clearly suggests a QP 
lifetime with the form l/r^p T. Finally, the peak in 
a I near Tc, most visible at low frequencies, has been at- 
tributed to phase fluctuations of the order parameter due 
to thermal vortex, anti- vortex pairs pl[ | . 

Following its success in YBCO, we attempt to describe 
the conductivity using a two fluid model, in which tr be- 
low Tc is comprised of condensate and QP contributions. 
The condensate's real, or dissipative, conductivity is a 
(^-function at = 0. Therefore, for any nonzero fre- 
quency CTi is due solely to the QPs. Following the anal- 
ysis of YBCO data, we choose a Drude conductivity for 
the QPs, aqp{uj,T)/c7Q = {pnTqp/h)/{l + uj^T^p). Here 
(Tq = / {fid) is the quantum conductivity of a stack of 
bilayers with spacing d = 15. 4A. The spectral weight, p„ 
is expressed in units of energy as will be described below. 

The two T dependent parameters, Tqp{T) and Pn{T) 
are clearly suggested by the data. Regarding Tqp (T) , the 
upper right panel of Fig. 1 suggests that I /Tqp ~ T over 
a wide temperature range. The form of Pn{T) follows 
from the T dependence of the superfluid density. Wc ob- 
tain ps (T) from the imaginary part of the conductivity, 
(72 , which is measured independently of cti in our experi- 
ment. We observe that Ps{T) = po — aT below about 30 
K, suggesting that p„ = aT at low T. For purposes of 
modeling, we have found that the most accurate descrip- 
tion of the terahertz conductivity results if we assume 
that pn — aT up to Tc. Of course the total normal fluid 
weight must increase faster than this in order that Pn~^Ps 
remains constant as T changes. However, we are mod- 
eling only the spectral weight contained in a Lorentzian 
component centered at zero frequency. 

In the right-hand panel of Fig. 1, we plot aqp{T), the 
Drude conductivity with Pn{T) and TqpiT) as described 
above. Focusing first on the highest frequencies in our 
range, we see that aqp (T) agrees rather well with the ob- 
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served ai{T). On the other hand, ai at low frequency 
is much larger than predicted by a Drude model for the 
QP conductivity. In fact, we have found it impossible 
to describe the conductivity in the range from 0.2-1 THz 
with a two-fluid model (p„(r) -I- Ps(T) is constant) re- 
gardless of the choice of Tqp{T). The difficulty is that 
the spectral weight in the low-frequency end of our range 
increases with decreasing T, rather than decreasing as 
expected for the normal fluid component in a two-fluid 
model. The conductivity spectra suggest the presence 
of an additional component, whose spectral weight in- 
creases with decreasing T. (Simulations based on p^ ] 
show that the spectral weight at low frequency cannot 
be caused by inadvertent coupling to the c-axis plasmon 
due to off-normal incidence). 
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FIG. 2. The spectral weight of the total conductivity from 
0.2-0.8 THz, E(r), the QP Drude conductivity, ^qp{T), and 
the difference between the two, S(T) — Eqp(T). The solid 
line shows the superfluid density, ps{T), multiplied by 0.18, 
showing the proportionality between this and the portion of 
the spectral weight not due to the QPs. 

To help identify the additional component of cr, we 
consider the difference between the measured spectral 
weight and the amount ascribed to QPs using the pa- 
rameters given above. Fig. 2 shows the spectral weight 
measured in our bandwidth, S(r), and that due to the 
normal fluid, Sgp(T), plotted versus T. Because of the 
extreme anisotropy of BSCCO we express spectral weight 
in terms of an areal density per Cu02 bilayer. The areal 
density is proportional to the phase or charge stiffness 
per bilayer, and can be conveniently expressed in ther- 
mal units of energy, or Kelvins. The spectral weights 
denoted by E are distinct from those discussed earlier 
labeled with p's. Whereas p{T) represents an integra- 
tion of a I over all frequencies, S(T) is confined to an 
integration over our experimental bandwidth. The dif- 
ference between S(T) and T,gp{T) is plotted as triangles. 
It is apparent that Sgp(T) is less than the total spectral 
weight observed. Notice that at low temperatures the dif- 



ference, S(T) — Eqp(T), is proportional to Ps{T), shown 
here multiplied by 0.18. The spectral weight in excess of 
T,qp is proportional to Ps only if we set l/Tgp{T) ~ T. 
If we choose l/Tgp{T) ~ with any (3 other than 1, 
the unaccounted for spectral weight does not have such 
a reasonable and recognizable T dependence. 

The proportionality of the residual spectral weight and 
Ps{T), suggests that the excess conductivity arises from 
fluctuations of the condensate order parameter. The fluc- 
tuations are not thermally generated because their spec- 
tral weight increases as T ^ 0. Although the temper- 
ature dependence is consistent with intrinsic quantum 
fluctuations, the associated conductivity is expected to 
be of order ctq or ^ 1.6 x 10^(ri~^ — m~^) which is 
far too small to explain the data. However, the fluctua- 
tion conductivity can be much larger in the presence of 
static (or quasistatic) spatial variation in the superfluid 
density. Two examples have been discussed recently in 
the literature. The first is a one-dimensional periodic 
modulation of the phase-stiffness parameter |Q. The 
second example is a "granular" superconductor in which 
the Josephson coupling between grains is randomly dis- 
tributed about its mean value 1 25 1 . Both of these models 
predict a contribution to ci above zero frequency which is 
not present if the phase stiffness is homogeneous through- 
out the material. Moreover, the spectral weight in this 
new contribution to ai is proportional to the mean value 
of the phase stiffness, or superfluid density, just as is ob- 
served in our experiments. In both cases the proportion- 
ality constant is {Aps)'^ / Ps'^ , the fractional mean square 
variation in p^. 

Motivated by these models, we attempt a decompo- 
sition of (Ti into QP and collective mode (CM) contri- 
butions, a I = Uqp + (Jem- The difference between cti 
and (jqp shows the collective mode contribution to be a 
low frequency peak contained, predominately, within our 
frequency range. We therefore test this decomposition 
with (Jem chosen to be a Lorenztian centered ai v = Q. 
This introduces two collective mode parameters: width, 
^cm{T), and spectral weight, Pcm{T). As is suggested by 
Fig. 2 we set the spectral weight of the collective mode, 
Pcm{T), equal to a fixed fraction, k, of ps{T). For the 
normal fluid spectral weight we choose p„ = aT as be- 
fore, and we vary Tcm and Tqp in order to find the best 
fit to the data. The best fit to l/Tqp{T) is shown in the 
upper right panel of Fig. 3. (The dashed line indicates 
the range where the large thermal fluctuation conductiv- 
ity dominates the QP contribution, making it difficult to 
extract Tqp). The width of the collective mode is found 
to be T independent within the noise level, with a value 
0.24 THz, and k = 0.30. 

The left panel of Fig. 3 compares cri(T) — (Jqp{T) 
with (Jcm{T), obtained using the parameter values given 
above. At each frequency (Ji{T) — aqp{T) is plotted as 
symbols and <Jcm{T) as a line. The difference between 
the total conductivity and that due to QPs is obviously 
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FIG. 3. Left panel: ai{T) - aqp{T) plotted at 0.2, 0.36 
and 0.64 THz as squares, circles and triangles, respectively. 
The lines show acm{T) from our model. Lower right panel: 
The difference between the data and the fit plotted on the 
same scale as the left panel, and ascribed to thermal phase 
fluctuations. Upper right panel: l/Tqp{T) for the QPs, due 
to this same model, plotted versus T with the T scale given 
by the panel below. 



decreasing temperature. We describe this contribution 
as a Lorentzian whose width is T independent and whose 
spectral weight is a constant fraction of PsiT). This ad- 
ditional dissipation could result from phase fluctuations 
in the presence of static or quasistatic spatial variations 
of the local superfluid density |2^,^^. From the spec- 
tral weight of this dissipation we infer a fractional mean 
square variation in ps of w 0.30. Adding this contribu- 
tion to that of QPs with 1 / r^p ^ T successfully describes 
the conductivity over the entire experimental range of 
frequency. Thus the transport lifetime in the supercon- 
ducting state increases with decreasing T far more slowly 
in BSCCO than in YBCO. The difference in lifetime may 
arise from the presence of same inhomogeneities which 
generate the collective mode dissipation. 
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and ONR Contract No. N00014-94-C-001. 



well described by the collective mode. What deviation 
there is between the model and the data is shown in the 
bottom right panel of Fig. 3. We see that difference 
corresponds to the phase fluctuations near ||2l| . 

At this point we comment on the collective mode con- 
ductivity spectrum. In both models described previously 
1^,^, the frequency of the collective mode is related 
to the screened plasma frequency of the condensate, or 
fu^p = pse^/ed. With e ^ 10, this yields fu^p ~ 1400 K, 
much higher than the energies probed in our experiment. 
However, the plasma mode can be strongly perturbed by 
thermal QPs. For example, the damping rate of the plas- 
mon, r, is given by fiT — PsOq /<yqp- For the values of Uqp 
seen in our experiment, this is of order 100 K, indicating 
that the plasma mode is highly overdamped. The corre- 
sponding conductivity will be centered at = and its 
width will depend on the damping and screening effects 
of the QP background. 

One final observation concerning our data is that 
causality dictates that the collective mode feature in 
(Ji{v) must also be clearly manifest in ij2{v). Since both 
(7i and (72 are measured directly and independently by 
our coherent spectroscopy this is an important check of 
our model's description of the conductivity. We find the 
total 172(1') to be a sum of contributions from the QPs, 
the condensate and the collective mode. <T2{v), therefore, 
is not described by a two-fluid model but rather by the 
model set forth above. 

In conclusion, we have shown that the dissipative con- 
ductivity of BSCCO in the superconducting state is not 
due solely to a normal fluid of QPs. There is an addi- 
tional contribution whose spectral weight increases with 
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